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Abstract 

We survey geometrical and especially combinatorial aspects of generalized Donaldson- 
Thomas invariants (also called BPS invariants) for toric Calabi-Yau manifolds, emphasizing 
the role of plane partitions and their generalizations in the recently proposed crystal melting 
model. We also comment on equivalence with a vicious walker model and the matrix model 
representation of the partition function Q. 



§ 1. Introduction 

The study of plane partitions, a three-dimensional generalization of partitions, 
has a long history of more than a century in mathematics [U [2]. There has recently 
been a renewed interest in this old topic, both among mathematicians and physicists 
alike, due to the pioneering discovery that topological A-model [3] on toric Calabi-Yau 
manifolds [I] can be described by a statistical mechanical model of plane partitions 
[51 [6] . In more mathematical language, plane partitions count Donaldson-Thomas (DT) 
invariants [8] , whose partition function is equivalent (9} [10] to that of the Gromov- 
Witten invariants under suitable parameter identifications. 

There is an interesting twist to this story, which is the topic of more recent stud- 
ies in this field. There we study "generalized Donaldson-Thomas invariants" [111 [12], 
which depend on moduli (mathematically stability conditions, or physically complexi- 
fied Kahler moduli El)- These "invariants" are invariant under a generic deformation of 
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moduli, but can jump when we cross real codimension loci (called walls of marginal sta- 
bility, which divide the moduli space into chambers) in the moduli space. This jumping 
is called the wall crossing phenomena, and there are general formulas |TTl [T2l [13] which 
govern this jumping phenomena. Generalized DT invariants are indeed generalizations 
of the original DT invariants, in the sense that the former coincide with the latter only 
in a specific chamber (hereafter called the topological string chamber) of the moduli 
space. 

Given the richness of these new invariants, the natural question is whether there 
are combinatorial counterparts to this geometric story. The goal of the present article 
is to provide an answer to this question. Generalized DT invariants on toric Calabi- 
Yau manifolds are described by a statistical mechanical model of "crystal melting" 
[T^t [T5| [16] , formulated here as an enumeration problem of plane partitions and their 
generalizations. We will also comment on the equivalence with a vicious walker model, 
following [17]. No prior knowledge in this field is assumed, and this paper is intended 
to be self-contained, at least as far as the combinatorial aspects are concerned. 

This article is organized as follows. In section 2 we define our combinatorial parti- 
tion function as a sum over suitable evolution of partitions. In section 3 we comment on 
the representation of the partition function as a unitary matrix integral. The derivation 
of this matrix integral is given in section 4, based on an equivalence with a vicious 
walker model. Appendix contains an introduction to generalized Donaldson-Thomas 
invariants, which readers can consult for geometric side of the story. 

§ 2. Definition of the Model 

We begin with the following theorem, which states that the partition function for 
generalized DT invariants Z g DT (see Appendix) on a toric Calabi-Yau manifold can be 
computed exactly by purely combinatorial methods: 

Theorem 2.1 (Szendroi [H], Mozgovoy-Reineke [15], Ooguri- Yamazaki |16j). 
For a toric Calabi-Yau manifold, the partition function for generalized DT invariants 
can be written as 

(2.1) -^gDT = -^crystab 

where the statistical mechanical model in the right hand side can equivalently formulated 
as (1) a crystal melting model (a configuration of molten atoms), (2) dimer model, or 
(3) a generalization of plane partitions (an evolution of partitions)^. 



The third formulation is available only when the toric Calabi-Yau manifold has no compact 4-cycles. 
All the examples in this paper fall into this category. 
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The goal of this section is to define Z crysta i, using the plane partitions and their 
generalizations. We use the third formulation, developed in [HI HH [201 ETJ • See [HI [15], 
[TBI [22] for the first and the second description. 

Remark. As we discussed in the introduction, the LHS of the equation (12. ID 
depends on the value of the moduli. Correspondingly, RHS also has moduli dependence, 
and we have different statistical mechanical models for different chambers of the moduli 
space. The schematic relation (12. ip should be interpreted this way. See (12.161) . 

Let us begin with standard notations. A partition A = (Aj) is a non- increasing 
sequence of integers Ai > A2 > . . . > such that A n = for sufficiently large n. The 
length |A| of a partition A is given by |A| := J2i -V We define a transpose A* by 

A*, :=#{j|A 2 < j}. 

This is simply a graphical transpose of a partition, as will be clear from the following 
example. 



n 



Example 2.2. For A = (4, 2, 1) = tr , |A| = 7 and A* = (3, 2, 1, 1) = □ 

Given two partitions A and //, we define A >- fi if and only if 

Ai = Hi + 1 or fXi for all i. 
We also denote A >~ \x if and only if A* >- /i*, or equivalently 

Ai > \i\ > A 2 > \iv > ■ ■ ■ ■ 



n 



Example 2.3. Two partitions A = (4, 2,1) = Lr and /u = (3, 2, 1) = 
both satisfy A >~ \i and A y 11. 

Now we define a plane partition (also called a 3d partition) as a sequence of parti- 
tions IT = {X(n)} ne z such that 

(2.2) ... 5 A(-2) 5 A(-l) 5 A(0) y A(l) y A(2) y . . . , 

and X(n) = {0} when \n\ sufficiently large. Define the length |n| of a plane partition 
IT = {A(n)} to be | XI | = J2 n l-^( n )l- Note that this is a finite sum by the assumption 
above. Our partition function is defined by 

(2.3) Z crystal (q) :=^ |n| - 

n 
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As is well-known, MacMahon's formula represents this as an infinite product 
(2.4) Z CTys Uq) = - q k y\ 

k>0 

which is the same as the generalized DT partition function for C 3 [3 J : El 

(2-5) ^crystal (<?)= ^Jdt (<?), 

For this reason we hereafter denote the LHS of the above equation (12. 5fl by tgl . 

Remark. In the definition of a plane partition ( 12. 2p we could have used y, -< 

instead of y, -<. This has the effect of replacing {X(n)} by {A(n)*}, and the partition 
function is the same either way. 

Remark. The choice of the weight in (1 2 . 3 j) is the same as in the Schur process of 



We have seen that plane partitions correspond to the simplest Calabi-Yau geometry 
C 3 . Our next task is to consider a set of partitions corresponding to the (resolved) 
conifold. We again consider a sequence of partitions IT = {A(n)}, but now with a 
slightly different interlacing conditions, with plus and minus appearing alternatingly: 

(2.6) . . . ^ A(-2) 5 A(-l) ^ A(0) y A(l) y A(2) y ... . 

For such a IT, we define |TT| := E n:e ven l^( n )l and similarly |II|i := Zln:odd l-M n )l- Tne 
conifold partition function is then defined by 

(0 7\ 7Conifold ( ' \ |n| |n|i 

\ Z ->) ^crystal l<?0, Ql) • — ^ % <?1 • 

II: satisfying (12.61) 



Then (12.11) in this example states that 

(2-8) ^t f a! d (90, qi) = Z^f old (q, Q) 

under the parameter identification 

(2.9) q = q qi, Q = q\. 

Infinite product expression for this partition function is known [3] from the study of the 
topological string: 

(2.10) Z™™ iold (q, Q) = M(q) 2 J](l + U (1 + q k Q~ l )\ 

k>0 k>0 



4 See Appendix for summary of these invariants. 

5 In the language of topological strings, q is related to the topological string coupling constant g s by 
q = —e~ 9s . 

6 In the topological strings, q = — e~ 9s , Q = — e~*, where g s is the topological string coupling 
constant, and t is the Kahler moduli of the resolved conifold. 
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The corresponding statement for ^crystal was shown combinatorially in |22j . 

The discussion above is a little bit imprecise because we did not specify the moduli 
dependence of the generalized DT invariants. We therefore consider the following more 
general partition function which includes the moduli dependence and applies to more 
general toric geometries. 

First, fix an integer L and a function p : {1/2, 3/2, . . . , L — 1/2} — > {±1}. We can 
periodically extend p to a map a : Z x / 2 — > {±1}, where 1>\/2 is a set of half-integers. 
This will determine a toric Calabi-Yau geometry. Second, to fix a moduli dependence, 
we define a bijection 9 : Z^ — > %i/2 such that0 

(2.11) 6(h + L) = 6(h) + L for all h e Z 1/2 , 

and 



2, 

i=i x ' i=i v 7 

We then define a generalized plane partition of type (L, p, 9) (whose totality we are 
going to denote by T>( l >p^) to be a sequence of partitions IT = {A(n)} such that 

ao0(i+l/2) ( 1\ 

A(i) < X(i + 1) for 6 I i + - J < 0, 
(2.13) ) 2 { 

A (i) >- A(z + 1) for ( i + - > 0. 



2, 

We define := J2 n =i mod L l^( n )l for z = 0, . . . , L — 1. We also define 



Qe- 1 (i-i/2)+i/2 • (i-i/2)+3/2 9e-i(i+i/2)-i/2 (0 1 (i-l/2)<0 1 (i + l/2)), 

^6>- 1 (i-l/2)-l/2 ' ?6)-i(i_i/2)-3/2 ^fl" 1 (t+l/2)+l/2 1 (i-l/2)>6 1 (z + l/2)), 



where we used qi for i G Z by extending 50, . . . , qL-i periodically, 

qi+L = qi- 

Note that q® =ld = qi. The partition function is defined by 

(2-14) Z^f2{qo,qi,...,q L -x)-.= £ (rf)' 11 ' ^)' 11 ' 1 • • • {A-i)^ 1 



7 This notation including half-integers looks cumbersome, but is useful to see the action of the Weyl 
group of the affine Kac-Moody algebra (281 [T9 , 20 . The parametrization by 6 actually covers only 
half of the chambers, but sufficient for our purposes here. The partition function becomes a finite 
product in the remaining half. 
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Example 2.4. Let us take L = I, p = +1 and then the only choice for 9 is 
9 = id, and dZH reduces to (Q. This means Z ( c ^ =+1 ' 9=id) = Z^ ystal . 

Example 2.5. Take L = 2, p(l/2) = +1 and p(3/2) = -1. can in general 
written as 

n n 1 1 33 

= 0„ : - t—T- ?i, — i — y — \- n. 

n 2 2 '2 2 

When 9 = 9q, (12 .131 ) is the same as (12. 6D . and the partition function (12 .14ft coincides 
with ( 12. 7| . The case of n 7^ corresponds to generalized DT invariants in other cham- 
bers. The corresponding partition function is given by [241 |25| 126] (again under the 
identification ( 12. 9| ) 

(2-15) ^t f ai ld (9o, ?r; n ) = M{qf J] (1 + g fc Q) fc J] (1 + q k Q~ 1 ) k . 

k>0 k>n 

In particular, in the limit n — > 00, this coincides with the commutative DT partition 
function for the conifold [Bj. In this limit our statistical mechanical model reduces to 
the gluing of two crystal corners (topological vertices) as in [H Hj . 

The general story goes as follows. We can construct from (L, p) a toric Calabi- 
Yau manifold X^ L,P \ which is one of the so-called generalized conifolds. This is a 
toric Calabi-Yau manifold without compact 4-cycles El, and has a connected string of 
L — 1 P s. Each P 1 is either a 0{— 2, 0)-curve or a 0(— 1, — l)-curve, depending on 
a(i - 1/2) = a(i + 1/2) or a(i - 1/2) = -a(i + 1/2). We can then consider the 
partition function of generalized DT invariants on X^ L,P \ 

The remaining task is to specify the moduli dependence, which in this case is given 
by an element of the Weyl group of the affine Lie algebra [28]. The corresponding 

partition function is denoted by Z^ T (q : Q; 9). Now the following theorem states that 
this partition function is the same as the crystal partition function of type (L, p, 9): 

Theorem 2.6 (Nagao [28]). 

(2-16) ^rystalteo, ■ ■ • , q L -l] 0) = Zg T (q, Q; 9), 

where the parameter identifications are given by 
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±<?o • • -<?l-i, Qi = ±qi (i = 1, . . • , L - 1). 



8 See [27] for discussion of Calabi-Yau geometries with a compact 4-cycle. 

9 This means that the overall sign change of p does not change the geometry. In (|2.13[) . this has the 

effect of replacing y,~<by>~,~<. As discussed previously in the case of C 3 , this does not change 
the partition function, but will change the matrix model representation of the partition function 
presented in the next section. 
10 The signs are determined from p. See [15] and section 3.5 of |29| . 
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Remark. We can consider further generalizations, by changing the boundary con- 
ditions at infinity. This generalized model counts "open generalized Donaldson-Thomas 
invariants" . See pH EDI EDI EE] • 

In the following we concentrate on the case of C 3 and the resolved conifold. 



§ 3. Matrix Model 

In the following sections we show that the crystal melting partition function defined 
in the previous section can be written as a unitary matrix integral: 

Theorem 3.1. 

(3.1) Z% UiGz)= lim / dUdetQ(U\q), 
where dU is the Haar measure of the unitary group and 

oo 

(3.2) e{u\q) = JJ(l+«g fc )(l + M"V +1 )- 

k=0 

Theorem 3.2 (Ooguri-Sulkowski-Yamazaki [1TJ ) . 

(3-3) ^crystal^QS Q't ®n) = ^n^matrbc^QS Q] 77,), 

where 

n 1 00 / 1 n~ 1 n k " 

^ c - - n (T^f n 

fc=l V y y fc=n+l V y 

and 

(3-5) Z™*£ d (q,Q;n) = lim / cft/det ( ^'f f[ (1 + Q^IT 1 ^) 

where the measure dU and the function Q(u\q) are the same as in Theorem 3.1. 

We shall give derivations of these results in the next section, but before going there 
some comments are in order. 

Remark. Theorem 13.11 and Theorem 13.21 for n = are not new, although n^O 
case has not previously appeared in the literature as far as the author is aware of. 
Theorem 3.1 seems to be well-known in the literature, and can be considered as a 
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reduction of a multi-matrix model in |32j. For Theorem 3.2 with n — 0, see the paper 
[33] , which also proves similar identities for other groups. A vicious walker model similar 
to the one for the conifold presented in the next section was also discussed in In 
their terminology >-,-< (>-,-<) are called "normal" ("super") time evolution, and the 
model is analyzed by identifies involving semi-standard Young tableaux and hook Schur 
functions (also called supersymmetric Schur functions). They also analyze the scaling 
limit of the model. See also |35|. 



Remark. The prefactor C n simplifies in the limit n = and n — > oo; Co = 1 and 
Coo = M(q). In particular in these cases C n is independent of Q. 

Remark. ^matrix, being a partition function of a unitary matrix integral, is a 
reduction of a r-function of a two-dimensional integrable Toda chain [36] (see e.g. [37] ) . 
Similar integrable structures appeared in topological strings context in [38]. See also 
[391 HO] for the appearance of thermodynamic Bethe Ansatz equations in the study 
generalized DT invariants. 

Finally, let us discuss the thermodynamic limit of our model, using the matrix 
integral given above. The thermodynamic limit is the limit g s — > 0, where the string 
coupling constant g s is related to the parameter q by q = e~ 9s 0- For small g s , the 
modular transformation of O with respect to g s gives 



e(e^|e- ffs ) =e~€l ■ (l + 0(e 



->3 



If we ignore non-perturbative terms in g s , this means that the matrix model ^atiiv 
reduces to a Hermitian Gaussian matrix model with unitary measure. This result was 
originally derived from the Chern-Simons theory on the conifold pTTj WI\ (see also 
for crystal melting description). 

The spectral curve for this Gaussian matrix model is given by the equation [12| | 

(3.6) e x + e v + e x - y ~ T + 1 = 0, 

where T = Ng s is the 't Hooft coupling. This is the mirror of the resolved conifold. In 
the limit of T — > oo (which we should take since — > oo), the curve reduces to 

(3.7) e x + e y + l = : 
which is the mirror of C 3 . 



Next we discuss the spectral curve for the conifold matrix model ( I3.5j) I I As before, 
we take the limit g s — > 0, N — > oo with T := Ng s fixed, but now we also take n — > oo, 



lx This is the parameter counting the size |I1| := y^^^ 1 of a generalized plane partition n. 
12 This is the spectral curve for the matrix model. Our statistical model can equivalently be written 
as a dimer model, which has its own version of the spectral curve. See [45j and [46J. 
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with r := ng s fixed. The spectral curve is given by [T7] 

(3.8) e x+y + e x + e y + Q ie 2x + Q 2 e 2y + Q 3 = 0, 



where 



Qi 



1 + fiQ 



(3.9) 



Q2 



(l + ^e 2 )(l + 2 )' 
1 + Qe 2 



Qs 



Q ■ 



(l + //Q)(l + //e 2 )' 
1 + 



(l + e 2 Q)(l + ^Q)' 



and \i = Q 1 q n ,e 2 = e T . This is the mirror [17] of the so-called closed vertex geometry, 
whose web diagram is shown in Figure [Tj. 



Figure 1. The web diagram for the close vertex geometry. Three IP's with size Qi,Q 2 , Q3 
appear symmetrically. 

There are two interesting observations on this result. First, ( 13. 9|) coincides with the 
mirror map for this geometry. Second, the curve ( I3.8[) is symmetric under exchanges 
of Q, [i = Q~ 1 q n and e 2 = e~ T . Namely, (1) the original Kahler moduli Q of the 
resolved conifold, (2) the chamber parameter n and (3) the 't Hooft parameter T ap- 
pear symmetrically in the spectral curve. This is an interesting result, which suggests 
a possible connection between continuum limit of the wall crossing formulas and the 
BCOV holomorphic anomaly equation |48j . 

In the matrix model we are interested in the limit of iV — > 00, which means T — > 00 
or equivalently e — > 0. With appropriate shifts of x and y, the equation (13. 8[) in this 
limit becomes 



The is the mirror of the so-called Suspended Pinched Point (SPP) geometry, with Q and 
\i being exponentials of flat coordinates representing sizes of its two P s, which encode 




(3.10) 



zy + e *+y + e * + (1 + Qfl ) e v + q = . 
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two copies of the initial 0(—l,—l) — >■ P 1 geometry, see figure [2] Not only does the 
spectral curve agree with the mirror curve of the SPP geometry in the limit of g s — > 0, 
but in fact the matrix integral reproduces the full topological string partition function 
all orders in g s expansion. Indeed, it is known that the SPP topological string partition 
function 



13 



(3.11) 



with Kahler parameters Q and /i, is equal to 

Qq k ) k (l- 



top 



H 

k=i y 



k\k 



_ g fc)3fc/2(i _ ^Qqk 



ifc" 



On the other hand, from the explicit structure of the BPS generating function and 
formulas ( 12. 15[) . ( I3.3[) . ( I3.4[) and ( I3.5p . we find that the value of the matrix integral, in 
the iV — > oo limit, is related to the above topological string partition function as 



(3.12) 



Zmatrix(<7, Q\ Tl) = Zf£? (q, Q,fJ,= Q 



-l qn] 



no 

k=i 



q 



k\k/2 



This result is consistent with the philosophy of the remodeling conjecture [39], which 
states that a set of invariants (symplectic invariants) [50] constructed recursively from 
the spectral curve coincide with topological string partition function on the same ge- 
ometry. Since symplectic invariants are defined by rewriting loop equations of matrix 
models purely in the language of spectral curves, the fact that the topological string par- 
tition function can be written as a matrix model would prove the remodeling conjecture. 
Indeed, this type of logic was used in [32], [51] to prove the remodeling conjecture for 
toric Calabi-Yau manifolds. It would be interesting to know whether similar recursion 
relations exist in other chambers. 





Figure 2. The web diagram for the SPP geometry. This geometry has two P^s with 
size Qi and Qi- 



§ 4. Derivation of the Matrix Models 

There are at least two derivations of the above-mentioned matrix models, one using 
the vertex operator formalism for free fermions [52] (see [T8l [T9l I2T] for discussion our 

13 In our context, a topological string partition function is a generalized DT partition function in 
the topological string chamber, where generalized DT invariants coincide with the original DT 
invariants of [7JIH]- This chamber is the analogue of 9 n=00 in the conifold 
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context) and another using the equivalence with a vicious walker model (non-intersecting 
paths). Both are presented in [17J. Here we comment on the latter method in the case 
6 = id. The derivation here is a slightly simplified version of the derivation in [17]. See 
also [32, 35, 51], which constructs similar matrix models in a particular chamber. In 
particular [51] treats arbitrary toric Calabi-Yau geometries. 

Let us fix a sufficiently large number N. This is the same iV for the size of the 
unitary matrix in the previous section, and in the end we take the limit iV — > oo. Define 

(4.1) h k {t):=\ N - k+1 (t) + k-l, (k = l,...,N). 
Since X(t) is a partition, we have 

(4.2) h k (t) < h k+ i(t), 
for all t. We also have the boundary condition, 

(4.3) h k (t) = k — 1 when \t\ large. 
Moreover, ( 12.21) means we have, for each step t, 

(4.4) h k (t + l)-h k (t) = 0oT -1, 
for t > and 

(4.5) h k (t + l) -h k (t) =0 or 1, 
for t < 0. 

Suppose that we fix a large positive (negative) integer t max (t m i n ). We are going to 
send these numbers of infinity. If we plot the value of {^fc(£)}t=t^ in for each k, we have 
a set of iV paths. Due to the conditions (j4.4|) . (14.51) the paths move on the graph shown 
in Figure [31 and ( 14. 3p means we have a fixed boundary condition. Finally, (14. 2 p means 
that iV paths are non-intersecting. Summing up, we have a statistical mechanical model 



of non-intersecting paths (also called a vicious walker model [S3] l 14 l). whose partition 
function is given by: 

(4.6) Z~ J2 l[q^ hi ^\ 

{hi(t)}: non-intersecting paths on the graph t 

where ~ shows that we neglected an overall multiplicative constant. 

At first sight it seems difficult in practice to implement the non-intersecting con- 
ditions for paths. The following theorem states that we can write the sum over non- 
intersecting paths as a determinant of a matrix, whose element is defined by a single 
path: 



14 We can also regards this model as a time evolution of N particles in one dimension. In this language 
the model is an exclusion process, a variant of the ASEP [54 1. 
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t=-3 t=-2 t=-1 t=0 t=1 1=2 t=3 



Figure 3. Top: an oriented graph for C 3 . Middle: an example of 3 non-intersecting 
paths shown as bold (red) arrows. The location of the k-th path at time t gives hk(t). 
Bottom: The corresponding evolution of Young diagrams. 
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Theorem 4.1 (Lindstrom [55], Gessel-Viennot [56]; Karlin-McGregor [57]). Sup- 
pose we are given an oriented graph without oriented loops. Suppose moreover that each 
edge e comes with a weight w(e). For a path p on the graph, we define w(p) to be the 
product of the weighs for all the edges on the path: w(p) := Y\ eep w{e) . We define F by 
summing over all non-intersecting paths {pi} (each pi starts from en and ends at bi): 

(4.7) F{{a i },{b i }) = 

{pi'.cii— >bi}: non-intersecting i 

and an N x N matrix G(ai,bj) by 

(4.8) G(a i ,b j ) = w{ 3>)- 

p: a path from ai to bj 

Then 

(4.9) F({a i },{b i })=det(G(a i ,b j )). 



Proof. When we expand the determinant det^j (G(a,i, bj)), we have contributions 
from non-intersecting as well as intersecting paths. However, contributions from the 
latter cancel out because they always come in pairs with an opposite sign (Figure 2]) 




Figure 4. When we expand det G, intersecting paths always come in pairs with an oppo- 
site sign. The reason is that we can exchange the label for paths after the intersection, 
without changing the paths themselves. 



Example 4.2. Consider an oriented graph with weights Wi, . . . ,Wq as shown in 
Figure EJ It is easy to see that 

(4.10) F({oi, a 2 }, {6i, b 2 }) = (wiw 4 )(w 3 w 5 ), 



15 When more than two paths intersect at a single point, we need to pick two of them according to a 
fixed ordering and apply the same argument. 
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and 

(4.11) G(ai,bj) = 

\W2W4 w 2 w G + w s w 5 
We indeed have detG = wiW4(w2Wq + W3W5) — (w2W4)(wiWq) = F. 

ai \n 1 ci w 4 bi 
•- 



a 2 w 3 c 2 w s b 2 

Figure 5. An example of an oriented graph. 

Remark. The fact that paths are non-intersecting is a manifestation of free 
fermions, and the determinant is interpreted as a Slater determinant. 

Therefore we have 

(4-12) Z^ ys Uq)=det(G l , J (q)) 1 

where Gij(q) is defined as a weighted sum over all possible paths which start at height 
i at t = £ m i n and end at height j at t = t max . As we can see from Figure El G(ai,bj) 
depend only on the difference i — j, and we thus have 

(4-13) Z^ ystal (q) = det(G l - J (q)). 

It turns out to be easier to write a generating function for G n than to write each 
G n separately: 

(4.14) f(z) := J2 = n^ 1 + II^ 1 + ^~ V +1 )- 



n 



Proof. To see this, note that a term in the expansion of the product is in one-to- 
one correspondence with a path. For example, for t < we take either 1 or zq l from the 
product, and the choice corresponds to the two possibilities in (14.51) . The change of the 
horizontal coordinates is measured by z, and taking the coefficient in front of z n means 
summing paths with height change n. The product in f )4.14p is over all non-negative 
integers n when we send t m j n — > — 00, t max — > 00. □ 



Finally, we have the following theorem: 
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Theorem 4.3 (Heine [58J, Szego [59]). Suppose thatf(z) — ^2 n G n z n . We then 

have 



(4.15) 




dU det f(U)= det Gj_,- 

1 <i,_7 < TV 



Remark. The RHS of the equation is often called a Toeplitz determinant of /. 

Proof. Diagonalize the unitary matrix U to be (e v/ ~ T< ^ 1 e^~^^ N ) . Then the 
integral J dU reduces to f Y[ dfa detij(e v/ ~ Tl ^ J ') deti > j(e~^~^ t ^), while the integrand 
becomes a product Yli /(e v/ ~ T ^ i ). After expanding the two determinants using the 
definition of the determinant, we can easily carry out the integral, and the result follows. 

□ 

This theorem, together with the form of f(z) in ( I4.14p . completes the derivation of 
the matrix model for C 3 . 



The analysis for the conifold is essentially the same, so let us summarize the result 
briefly. By defining hk(n) again as in (j4.ll) . we again have (14.21) and (14.31) . except that 
(14.41) . ( 14. 5ft are going to be replaced by 



1. When t is odd, 



(4.16) h k (t + l)-hk{t) 



0, 1 (t<0), 
0,-1 (t > 0). 



2. When t is even, 

(4.17) ... < h k -!(t + l) < h k (t) < h k (t + l) < h k+1 (t) < .... 
for t < and 

(4.18) . . . < fc fc _i(t) < h k (t + 1) < h k {t) < h k+1 (t + 1) < . . . . 
for t > 0. 

These conditions mean {h k (n)} move on the graph shown in Figure [61 Note that the 
structure of the graph is different depending on whether t is even or odd. 
Again by using Theorems 14.11 and 14.31 we have 

(4.19) Z = [ dU det f(U), 



16 



Masahito Yamazaki 




t=-3 t=-2 t='-1 t=0 t=1 t=2 t=3 





V. 










,X 










Si 


X 






V. 




X 


X 








A 




ft- 


,N 


>■ 






V, 






X 



t='-3 t=-2 t="-1 t=0 t=1 t=2 t=3 




Figure 6. Top: an oriented graph for the conifold. Middle: an example of 3 non- 
intersecting paths on the graph shown in red. Bottom: the corresponding evolution of 
Young diagrams. 
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where 

(420 f(u)=r[ i+iMigf n 1+ ^ n+lz ~ 1 

This is nothing but the expression ( 13.5 j) for n = 0. 
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§ Appendix A. Generalized Donaldson-Thomas Invariants 

In this appendix we briefly summarize the ingredients of the generalized Donaldson- 
Thomas invariants. Our discussion in this section is far from rigorous and at best 
schematic, since the main focus of this paper is more on combinatorial aspects presented 
in the main text. 

For the definition of generalized DT invariants, we need the following: 

• X: a Calabi-Yau 3- fold. 

• a "charge lattice" : 

H even (X; Z) = H (X; Z) © H 2 (X; Z) © H 4 (X; Z) © H 6 (X; Z). 

• complexified Kahler moduli of X: 

U = Bi + V^lh, (i = l,...,dim# 2 (X;Z)), 

where the real part Bi (imaginary part fcj) denotes the B-field flux through (the 
volume of) the i-th 2-cycle. 

• A central charge function Z 1 (t), which depend on t := {ti} and linearly on 7 G 
H even (X; Z) |_. 

16 This is part of the data for the stability conditions [60J. 
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With these data we can define "generalized DT invariants" 

0( 7 ; t) e Q. 

For concreteness, in this paper we restrict ourselves to the following situations: 

• X: a toric Calabi-Yau 3-fold without compact 4-cycles. For example, X can be 
C 3 or the (resolved) conifold. 
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• In the charge lattice H even (X; Z) we only consider the following set of charges 

H even (X; Z) = H (X; Z) © H 2 (X; Z) © H 4 (X; Z) © H 6 (X; Z) 
7 =( n, 13 = {Pi}, 0, 1 ). 

We then have a set of integer invariants 

O( 7 = (n,/3,o,i);t) ez. 

Instead of studying these invariants separately, it is useful to define their generating 
function: 

(Appendix A.l) Z gDT (g, Q- 1) = ^ Q (7 = (n, /?, 0, 1); t) q n Q^ 

n,/9 

where Q := {Qi} denotes a set of parameters and Q@ := YliQi - This is the partition 
function for generalized DT invariants studied in the main text. 
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